Shear-rate dependent transport coefficients for inelastic Maxwell models 



Vicente GarzcB 

Departamento de Ftsica, Universidad de Extremadura, E-06071 Badajoz, Spain 

(Dated: February 1, 2008) 

The Boltzmann equation for d-dimensional inelastic Maxwell models is considered to analyze 
transport properties in spatially inhomogeneous states close to the simple shear flow. A normal 
solution is obtained via a Chapman-Enskog-like expansion around a local shear flow distribution 
j' ' that retains all the hydrodynamic orders in the shear rate. The constitutive equations for 
the heat and momentum fluxes are obtained to first order in the deviations of the hydrodynamic 
field gradients from their values in the reference state and the corresponding generalized transport 
coefficients are exactly determined in terms of the coefficient of restitution a and the shear rate a. 
Since f (0 ^ applies for arbitrary values of the shear rate and is not restricted to weak dissipation, the 
transport coefficients turn out to be nonlinear functions of both parameters a and a. A comparison 
with previous results obtained for inelastic hard spheres from a kinetic model of the Boltzmann 
equation is also carried out. 
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I. INTRODUCTION 



Granular gases are usually modelled as a gas of hard spheres dissipating part of their kinetic energy during collisions. 
In the simplest model, the grains are taken to be smooth so that the inelasticity of collisions is characterized only 
, t . through a constant coefficient of normal restitution a < 1. For a low-density gas, the Boltzmann kinetic equation has 
been conveniently modified to account for inelastic binary collisions p], Q and the Navier-Stokes transport coefficients 
H for states with small hydrodynamic gradients have been computed by means of the Chapman-Enskog method Q 
around the local version of the homogeneous cooling state. As in the case of elastic collisions, the transport coefficients 
are given in terms of the solutions of exact linear integral equations which are approximately solved by considering the 
leading terms in a Sonine polynomial expansion of the velocity distribution function. On the other hand, in spite of 
this approach, the Sonine predictions compare in general quite well with computer simulations, even for quite strong 
dissipation conditions [f|. 

Needless to say, the mathematical difficulties of solving the Boltzmann equation for inelastic hard spheres (IHS) 
increase considerably when one considers situations for which large gradients occur and more complex constitutive 
^vq ' equations than the Navier-Stokes ones (which are linear in the spatial gradients) are required. For this kind of 
\ situations one has to resort to alternative approaches, such as the use of simplified kinetic models, where the true 
Boltzmann collision operator is replaced by a simpler collision model (e.g., the BGK model [6]) that preserves its 
most relevant physical features. In the case of elastic collisions, this route has been shown to be very fruitful since 
ly-j , most of the exact solutions derived from the BGK model agree rather well 0] with results obtained by numerically 
solving the Boltzmann equation by means of the direct simulation Monte Carlo (DSMC) method Q. Nevertheless, 
much less is known for inelastic collisions although some exact results obtained in the simple or uniform shear flow 
(USF) problem 0] and in the nonlinear Couette flow state [l(| support the reliability of the inelastic version of the 
BGK model @ for granular gases as well. The USF state is perhaps the simplest flow problem since the only nonzero 
hydrodynamic gradient is du x /dy = a = const, where u is the flow velocity and a is the constant shear rate. Due 
to its simplicity, this state has been widely studied in the past for elastic and inelastic gases 11 1 as an ideal 



testing ground to shed light on the intricacies associated with the response of the system to strong shear rates. Very 
recently [lH, the inelastic BGK model has been applied to determine the transport coefficients for a granular gas in 
spatially-inhomogeneous states close to the USF. The heat and momentum fluxes were evaluated to first order in the 
deviations of the hydrodynamic field gradients from their values in the reference USF state. Given that the system 
is strongly sheared, the corresponding transport coefficients turn out to be nonlinear functions of both the shear 
rate and the coefficient of restitution. This is the main new ingredient of these constitutive equations. In addition, 
due to the mathematical difficulties involved in the general problem, results were restricted to a particular kind of 
perturbations for which the steady state conditions of the USF hold [l2| ■ This allowed us to perform a linear stability 
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analysis of the hydro-dynamic equations with respect to the USF state. 

As happens for elastic collisions, a possibility to analyze transport around USF retaining the explicit form of the 
Boltzmann collision operator is to consider inelastic Maxwell models (IMM), i.e., models for which the collision rate 
is independent of the relative velocity of the two colliding particles. Although these IMM do not correspond to any 
microscopic interaction, it has been shown by several authors [131 ] that the cost of sacrificing physical realism can be 
in part compensated by the amount of exact analytical results derived from this interaction model. This is the main 
reason why IMM have attracted the attention of physicists and mathematicians since the beginning of the century. 
On the other hand, beyond its academic interest, it must be remarked that recent experiments [l4j for magnetic grains 
with dipolar interactions in air in a two-dimensional geometry have been accurately described by IMM. The excellent 
quantitative agreement found between the magnetic particles experiments [l4| and IMM suggests that magnetic 
particles can be considered as an ideal experimental probe for the predictions of this analytically tractable kinetic 
theory. 

The main advantage of using IMM instead of IHS is that a velocity moment of order k of the Boltzmann collision 
operator only involves moments of order less than or equal to k. This allows one to evaluate the Boltzmann collision 
moments without the explicit knowledge of the velocity distribution function [Til ] . Thanks to this simplification, the 
velocity moments of IMM under USF can be computed as functions of the shear rate. Very recently [16[ , the first 
nontrivial shear-rate dependent moments of USF for IMM have been explicitly computed: the second-degree moments 
which are directly related to rheological properties and the fourth-degree moments whose knowledge is necessary to 
get the heat flux around USF. The knowledge of these moments allows one to reexamine the problem studied in 



Ref. [12[ in the context of the Boltzmann equation and without any restriction on the kind of perturbations to USF 



considered. In addition, the comparison between the transport coefficients derived here with those obtained from IHS 
[l2| can be used again as a test to assess the degree of usefulness of IMM as a prototype model for granular flows. 
Previous comparisons between IMM and IHS for inhomogeneous situations [Tt], [H, Gil have shown in general good 
agreement (especially for low order moments), confirming the reliability of IMM as a good approximation to model 
granular systems. 

Since I am interested in studying heat and momentum transport in a strongly sheared granular gas, the physical 
situation is such that the gas is in a state that deviates from the USF by small spatial gradients. Under these 
conditions and taking the USF state as the reference one, the Boltzmann equation is solved by means of a 
Chapman-Enskog-like expansion [2(| around the distribution f(°\ Since the latter applies for arbitrary values of the 
shear rate a, the successive approximations in the Chapman-Enskog expansion will retain all the hydrodynamic orders 
in a. Therefore, the non-equilibrium problem studied here deals with two kind of spatial gradients: small gradients 
due to perturbations of the USF and arbitrary large gradients due to the background shear flow. In this paper, the 
calculations will be restricted to first order in the spatial gradients of density, temperature and flow velocity. At this 
level of approximation, the momentum transport is characterized by a viscosity tensor r]ijk£, while the heat flux is 
expressed in terms of a thermal conductivity tensor Kij and a new tensor /z^- (not present in the elastic case). The 
set of generalized transport coefficients rjijkti and fJ-ij are nonlinear functions of a and a. The determination of 
this dependence for IMM is the primary target of this paper. 

The plan of the paper is as follows. In Sec. |TT1 the Boltzmann equation for IMM is introduced and a brief summary 
of relevant results derived for the USF problem is given. The Chapman-Enskog expansion around USF is described 
in Sec. IIIII and the linear integral equations defining the generalized transport coefficients are displayed. Section IIVI 
deals with the explicit evaluation of the transport coefficients associated with the momentum and heat fluxes. The 
details of the calculations are displayed along several Appendices. The dependence of some of the above coefficients 
on the shear rate and on the coefficient of restitution is illustrated and compared with known results obtained for IHS 
[T2I ]. The comparison shows in general qualitative good agreement, especially in the case of the viscosity tensor rjijke- 
The paper is closed in Sec. IVII with a brief discussion on the results reported in this paper. 



II. INELASTIC MAXWELL MODELS AND UNIFORM SHEAR FLOW 



Let us consider a granular fluid modelled as an inelastic Maxwell gas. The inelasticity of collisions among all pairs 
is accounted for by a constant coefficient of restitution < a < 1 that only affects the translational degrees of freedom 
of the grains. At a kinetic level, all the relevant information on the state of the system is provided by the one-particle 
velocity distribution function /(r, v, t). In the low density regime the inelastic Boltzmann equation [HQ gives the time 
evolution of /(r, v,i). The corresponding Boltzmann equation for inelastic Maxwell models (IMM) can be obtained 
from the inelastic Boltzmann equation for inelastic hard spheres (IHS) by replacing the rate for collisions between 
two particles (which is proportional to the relative velocity in the case of IHS) by an average velocity-independent 
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collision rate. With this simplification and in the absence of an external force, the Boltzmann equation for IMM reads 

d 



H +vV\f(r,v,t) = J[v|/(i), /(*)], (2.1) 

where the Boltzmann collision operator is 

J [ v i|/'/] = ^/ dv 2 / ^[""'/(viJ/WJ-ZlviJ/W] • (2.2) 

Here, n is the number density, = 2n d / 2 /T(d/2) is the total solid angle in d dimensions, <x is a unit vector along 
the line of the two colliding spheres, and g = vi — V2 is the relative velocity of the colliding pair. In addition, the 
primes on the velocities denote the initial values {v'^Vj} that lead to {vi,V2} following a binary collision: 

vi =v 1 -i(l + a- 1 )(?-g)?, v 2 =v 2 + i(l + a- 1 )(?-g)& . (2.3) 

The effective collision frequency u) is independent of velocity but depends on space and time through its dependence 
on density and temperature. Here, I will assume that to oc nT q , with q > 0. The case q = will be referred here 
as Model A, while the case q ^ will be called Model B. Model A is closer to the original Maxwell model of elastic 
particles, while model B is closer to hard spheres when q = i. Furthermore, the collision frequency u> can be also seen 
as a free parameter of the model, determined to optimize the agreement with some property of interest of the original 
Boltzmann equation for IHS. As in previous works [l2i [II, Gjj] , u> is chosen here to guarantee that the cooling rate of 
IMM be the same as that of IHS (evaluated at the local equilibrium approximation) [ItJ • With this choice, one gets 

w=^y^b, v =A(q)nT«, (2.4) 

where the value of the quantity A{q) is irrelevant for our purposes. Henceforth, I will take this choice for u>. 

There is another more refined version of IMM [2l[ where the collision rate has the same dependence on the scalar 
product (<r • g) as in the case of IHS. However, both versions of IMM lead to similar results in problems as delicate 
as the high energy tails and so I will consider here the simplest version given by Eqs. (|2.2p and (|2.3[) . 

The first d + 2 velocity moments of / define the number density 



the flow velocity 



n(r,t) = / dv/(r,v,t), (2.5) 



(2.6) 



n{r,t) J 



and the granular temperature 

T M) = ^b)/ ^ 2 M)/(r,v,i), (2.7) 

where V(r,i) = v — u(r, t) is the peculiar velocity. The macroscopic balance equations for density, momentum, and 
energy follow directly from Eq. (|2.ip by multiplying with 1, rav, and ^mv 2 and integrating over v: 

D t n + nV -u = , (2.8) 
D t Ui + (mn)~ l VjPij = Q , (2.9) 



D t T + 1. (V ■ q + PijVjUi) = -(T , (2.10) 

dn 



where D t = dt + u • V. The microscopic expressions for the pressure tensor P, the heat flux q, and the cooling rate ( 
are given, respectively, by 



P(r,t)= / dvmWf(r,v,t), (2.11) 
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q(r, t) = J rfv ±mV 2 V /(r, v, t), (2.12) 

C(M) = -^^/^J[r,v|/(t)]. (2.13) 

The balance equations (|2.8p - (|2.10p apply regardless of the details of the interaction model considered. The influence 
of the collision model appears through the a-dependence of the cooling rate and of the momentum and heat fluxes. 
In particular, the cooling rate £ is given by [ItJ 

1 - a 2 Gf + 2., 9 , , n<1 , 
C=^ r « = - £r (l-a?V (2.14) 

Let us assume that the gas is under USF. This idealized macroscopic state is characterized by a constant density, 
a uniform temperature, and a simple shear with the local velocity field given by 

Ui = aijTj, ciij = a5 lx S jy , (2-15) 

where a is the constant shear rate. This linear velocity profile assumes no boundary layer near the walls and is 
generated by the Lees-Edwards boundary conditions (24j | , which are simply periodic boundary conditions in the local 
Lagrangian frame moving with the flow velocity (25|. Since the heat flux is zero in the USF problem, the balance 
equation for the energy (|2.10[) reads 

v -^d t \nT = -Q* ~^P* yl (2.16) 

where = C/^Oi a * — a /^0i P£ y — Pxy/p, P = nT being the hydrostatic pressure. Equation (|2. 16[) shows that the 
temperature changes in time due to the competition between two (opposite) mechanisms: on the one hand, viscous 
(shear) heating and, on the other hand, energy dissipation in collisions. The reduced shear rate a* is the nonequilibrium 
relevant parameter of the USF problem since it measures the departure of the system from equilibrium. Note that, 
except for Model A (q — 0), a*(t) oc T(t)~ q is a function of time through its dependence on temperature. Since in 
the hydrodynamic regime P xy {t) depends on time only through its dependence on a*(t) then for q ^ a steady 
state is eventually reached in the long time limit when both viscous heating and collisional cooling cancel each other 
and the fluid autonomously seeks the temperature at which the above balance occurs. In this situation, a* and a are 
not independent quantities but they are related through the steady state condition: 



a*P: v = ~C- V2.X7) 



However, when q = 0, the collision frequency i/q is independent of temperature and a* remains constant in time, so 
that there is no steady state (except if a* takes the specific value given by (|2.17[) ). Consequently, only in the case of 
Model A the reduced shear rate a* and the coefficient of restitution a are in general independent parameters in the 
USF state. Note that, although the temperature changes in time, the distribution of velocities relative to the thermal 
speed ^/Wj m is expected to reach, after a sufficient number of collisions per particle, a stationary form that only 
depends on the control parameters a* and a [22| . As a consequence, when the velocity moments are conveniently 
scaled with the thermal speed, they reach stationary values after a kinetic transient regime [H, [23| . These steady 
values are nonlinear functions of both the (reduced) shear rate and the coefficient of restitution. More details on the 
USF state for dissipative systems can be found in Ref. [l6j . 

The USF problem is perhaps the nonequilibrium state most widely studied in the past few years both for granular 
and conventional gases 0, [H| • At a microscopic level, it becomes spatially homogeneous when the velocities of the 
particles are referred to the Lagrangian frame of reference co-moving with the flow velocity u [25]. In this frame, the 
one-particle distribution function adopts the uniform form, /(r, v) — » /(V), and the Boltzmann equation (|2.ip reads 



dt-aV y ^ r )f(V) = J[V\fJ] . (2. IS! 



Upon writing Eq. (|2.18p use has been made of the identity 



(2.19) 
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where in the last step I have taken into account that / depends on u through the peculiar velocity V. The elements of 
the pressure tensor provide information on the rheological properties of the system. These elements can be obtained 
by multiplying the Boltzmann equation (|2.18p by mViVj and integrating over V. The result is (TH IT5. 

d t P i3 + auPjt + ajiPu = -v Q \ 2 {Pij - pSij) - (pSij, (2.20) 

where 

^o|2=C + ^±^o. (2.21) 

In the case of Model A (q = 0), the set of first-order differential equations (12.201) can be exactly solved [HI]. In 
terms of the reduced elements P*a = Pij(t)/p(t), the solution can be written as 

_ 1 + 2d 7 (a) _ _ 1 (9 99] 

i + 2 7 (s) ' yy zz i + 27(a) ' [ ' 



where 



P* _ rf 7(a) _ a 
xy a [1 + 2 7 (S)] 2 



a . (l + «) 2 
a= , ^o|2 = ^o|2-C = ; v Q , 

W |2 4 



and 7(a) is the real root of the cubic equation 



(2.23) 



(2.24) 



7 (1 + 2 7 ) 2 



namely, 



7(3) = 3 sinh2 



- cosh Ml + —a 2 
6 \ a 



Insertion of Eqs. (|2T2^|l and fFZ3\i into Eq. (|2~2T^) yields 

9 t lnT = 2A, 

where I have called 

A = 27(a)w |2 - C- 



(2.25) 



(2.26) 



(2.27) 



(2.28) 



In the elastic limit (a = 1), a = a*, and Eqs. (|2.22p - (|2.28l) reduce to the well-known exact solution obtained long time 
ago by Ikenberry and Truesdell for Maxwell molecules [2Q|. Equation (|2.27[) shows that T(t) either grows or decays 
exponentially. The first situation happens if 27(a)a>o|2 > C- I n that case, the imposed shear rate is sufficiently large 
(or the dissipation in collisions is sufficiently low) as to make the viscous heating effect dominate over the inelastic 
cooling effect. The opposite happens if 27(a)w |2 < C- 

Beyond rheological properties, the next nontrivial moments in the USF problem are the fourth-degree velocity 
moments. These moments are needed to determine the transport properties of the gas in states close to USF, which 
is the main goal of this paper. Very recently, the dependence of the fourth-degree moments on the shear rate and 
the coefficient of restitution has been explicitly obtained for IMM [l6(. Their explicit expressions are displayed in 
Appendix [X] for a three-dimensional gas (d = 3) . In a similar way to the case of elastic Maxwell molecules Q , it has 
been shown that, for a given value of a, the fourth-degree moments are divergent for shear rates larger than a certain 
critical value a c (a). This singular behavior of the moments reflects the existence of an algebraic high- velocity tail in 
the distribution function. However, for practical reasons, since in general the numerical value a c (a) is rather large, 
nonlinear shearing effects are still significant for a < a c . 

As said before, in the case of Model B (q ^ 0), after a transient regime the system achieves a steady state, so that 
A = and Eq. (pT2"8) leads to 



7(a) 



a 



d + 21 



2d 1 + a' 



(2.29) 
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The steady solution for the pressure tensor is still given by Eqs. (|2.22[) and (|2.23[) , except that a (or, equivalently a*) 
and a are not independent quantities. The explicit dependence of the steady-state value a* (a) of the reduced shear 
rate can be easily obtained by inserting the condition (|2.29|) into Eq. (|2.25|) : 

*/ \ /rf + 2 ~d+l — a , n 

= y^-(l-a 2 ) — . (2.30) 

It must be remarked that the steady-state results are "universal" in the sense that they hold both for Model A and 
Model B, regardless of the precise dependence vq(£) For elastic collisions, Eq. (|2.30[) yields a* — and so the 

equilibrium results are recovered, i.e., P*j = 5{j. The analytical results obtained for the steady USF state in the case 
of Model B [l8[ agree quite well with Monte Carlo simulations of the Boltzmann equation for IHS [13, HH , even for 
strong dissipation. 



III. SMALL PERTURBATIONS AROUND THE UNIFORM SHEAR FLOW STATE 

Let us assume now that we disturb the USF by small spatial perturbations. The response of the system to these 
perturbations gives rise to additional contributions to the momentum and heat fluxes, which can be characterized 
by generalized transport coefficients. Since the system is strongly sheared, the corresponding transport coefficients 
are highly nonlinear functions of the shear rate. The goal here is to determine the shear-rate dependence of these 
coefficients for IMM. 

To analyze this problem, one has to start from the Boltzmann equation (|2.ip with a general time and space 
dependence. First, it is convenient to keep using the relative velocity V = v — uo, where uo = a ■ r is the flow velocity 
of the undisturbed USF state. Here, the only nonzero element of the tensor a is a^- = a5i X 5j y . On the other hand, in 
the disturbed state the true velocity u is in general different from Uo, i.e., u — Uq + Su, 5u being a small perturbation 
to uo. As a consequence, the true peculiar velocity is now c = v u = V <5u. In the Lagrangian frame moving with 
velocity u , the convective operator v ■ V can be written as 



v • V/ = (V + uo) • V/ = -aV y ^f + (V + u ) • V/, (8.1 ) 



where the derivative V/ in the last term must be taken now at constant V. According to the identity (|3.ip . the 
Boltzmann equation reads 



9 f T/ 9 

dt f - aV *dv> 



f - aV v7^rf + (V + u ) • V/ = J\V\f, /]. (3.2) 



I am interested in computing the transport coefficients in a state that slightly deviates from the USF. For this reason, 
I assume that the spatial gradients of the hydrodynamic fields 

A(r, t) = {n(r, t), T(r, t), 6u(r, t)} (3.3) 

are small. Under these conditions, a solution to the Boltzmann equation (|3.2p can be obtained by means of a 
generalization of the conventional Chapman-Enskog method [3| , where the velocity distribution function is expanded 
about a local shear flow reference state in terms of the small spatial gradients of the hydrodynamic fields relative to 
those of USF. This type of Chapman-Enskog-like expansion has been considered in the case of elastic gases to get the 
set of shear-rate dependent transport coefficients @, in a thermostatted shear flow problem and it has also been 
recently considered [12], |2fJ, |3fJ in the context of inelastic gases. 

The Chapman-Enskog method assumes the existence of a normal solution in which all space and time dependence 
of the distribution function occurs through a functional dependence on the fields A(r,t), i.e., 

/(r,V,i) = /[VL4(r,t)]. (3.4) 

This solution expresses the fact that the space dependence of the shear flow is completely absorbed in the relative 
velocity V and all other space and time dependence occurs entirely through a functional dependence on the fields 
A(r, t). This functional dependence can be made local by an expansion of the distribution function in powers of the 
hydrodynamic gradients: 

/(r, V,i) = fM(V\A(r,t)) + f^(V\A(r,t)) + ■■■, (3.5) 

where the reference zeroth-order distribution function corresponds to the USF distribution function but taking into 
account the local dependence of the density and temperature and the change V — > c = V — Su(r, t). The new feature 
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of this Chapman-Enskog expansion (in contrast to the conventional one) is that the successive approximations 
are of order k in the gradients of n, T, and Su but retain all the orders in the shear rate a. More technical details on 
this Chapman-Enskog- like type of expansion can be found in Appendix [5] 

The expansion (|3.5[) yields the corresponding expansion for the fluxes and the cooling rate when one substitutes 
([33]) into their definitions (|23Tj) - (|2~T3)) : 



p = P (0) + p (D + ... ! q = q (0) +q (D + ... i C = C (0) +C (D + ... _ (36) 
Finally, as in the usual Chapman-Enskog method, the time derivative is also expanded as 

d t = df ] + d { t 1] + d[ 2) + • ■ • , (3.7) 

where the action of each operator is also defined in Appendix |B] In this paper, only the zeroth and first order 
approximations will be considered. 

A. Zeroth-order approximation 

Substituting the expansions (|3.5p and (|3.7p into Eq. (|3.2p . the kinetic equation for is given by 

9 t (0) / (0) -^^ r /(°)= J[V|/(°),/(°]. (3.8) 



To lowest order in the expansion the conservation laws give 

~ dn^ xy 



d^n^O, 9 f) T= _ ^ p(0)_ rc (3 _ 9) 



df>8v,i + aijSuj = 0. (3.10) 

Upon writing the second identity in Eq. (|3.9p I have taken into account that the effective collision frequency to cx nT q 
is assumed to be a functional of / only through the density and temperature. Consequently, u^ ' = u>, lu^ ~ lu^ — 
■ ■ ■ = and, using Eq. (|2T4|) . C (0) = (, C (1) = C (2) = • • • = 0. It must be noticed that, in the case of IHS, is 
different from zero but quite small [12j. 

Since is a normal solution, its dependence on time only occurs through n, Su and T: 

Q(0M0) 9 / (0) fl(0) 3/<») („, ()/(») (0). 



where I have taken into account again that f ( °^ depends on <5u through c. Substitution of Eq. (|3.11| into Eq. ([3 
yields the following kinetic equation for /'°); 

- (1-aPM +T C ) A/(o) _ ac,A /( o) = J[V | / (o) ) /(P]> (3 . 12) 

The zeroth-order solution leads to q^°) = 0. In the case of Model A, the zeroth-order pressure tensor is given by Eqs. 
(12.221) and (|2.23|) while the expressions of the fourth-degree moments are displayed in Appendix^] As shown in Refs. 
[l2Ll20L[30l|. for given values of a and a, the steady state condition (|2.17ll establishes a mapping between the density 
and temperature so that every density corresponds to one and only one temperature. Since the density n(r, t) and 
temperature T(r, t) are specified separately in the local USF state, the viscous heating only partially compensates 
for the collisional cooling and so, <9 t T ^ [20|. Consequently, the zeroth-order distribution /(°) depends on time 
through its dependence on temperature and the quantities a* and a must be considered as independent parameters 
for general infinitesimal perturbations around the USF state. This fact gives rise to new and conceptual practical 
difficulties not present in the previous analysis made for elastic thermostatted gases [29| . 
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The set of equations for P^ follows from Eq. (|3.12|) as 



—aP^+Tc) — P (0) 
dn xy + ^ dT » 



auP, 



(0) 



ajiP t 



(0) 



v , fP (0) 
- v 0\2\^ij 



p5ij) - C,p8i 



(3.13) 



The dependence of P^ on temperature occurs explicitly through the hydrostatic pressure p — nT and through its 
dependence on a*. Consequently, 



d 



d 



dT 



dT 1 11 



d 



da" 



where here P* = P^ /p. Thus, in dimensionless form Eq. (|3.13|) becomes 



^ P xy + C 



1-qa*^- I P t 
da* ' 



an Pa. 



AM 



(3.14) 



(3.15) 



where vt 2 = v a \ 2 jv a . For q = (Model A), the solution to Eq. (|3~l"5|) is given by Eqs. (f2T22Tl and (f2~2"3"f . For 
g ^ (Model B), Eq. (|3.15[) must be solved numerically to get the dependence of P*= on a*. In the case of IHS 
(q = tj), a detailed study of the dependence of P*j on a* has been carried out in Ref. [23|]. For elastic gases (a = 1), 
a comparison between the results derived for Maxwell molecules (q = 0) and hard spheres (q = h) at the level of 
rhcological properties has shown that both results are very close [7J so that there is a weak influence of the interaction 
model on transport properties. This suggests to expand the transport properties in powers of the interaction parameter 
q as an alternative to get accurate analytical results for non-Maxwell molecules. For inelastic collisions, I expect that 
such a good agreement is also kept and so the following formal expansion in q is considered: 



P* j (a*) = P* j , A (a*)+q&P* j (a*) + 



(3.16) 



where P*j ^(a*) is the known result for Model A. Inserting Eq. (|3.16[) into Eq. (|3.15[) and neglecting terms nonlinear 
in g, one gets 



AP* 



6A* 



7 



^ |2 (l + 2 7 )(l + 6 7 ) 2 



(3.17) 



X*a 1 - 67 

*y~ ^ |2 (l + 2 7 )2(l + 6 7 ) 2: 



(3.18) 



where co^ 2 — wq^/^o an d A* = V^o, A being given by Eq. (|2.28|) . Upon writing Eqs. (|3.17[) and (|3.18[) use has been 
made of the relations 







P* — 

"„„ A — 



4", 



da* VV ' A (I + 27XI + 67)' 
a da* xy ' A a (l + 2 7 ) 2 (l + 6 7 )' 



(3.19) 



(3.20) 



From Eq. (|3.16|) and taking into account Eqs. (|2.22p and (|2.23p . the expressions for the elements of the pressure 
tensor for Model B can be written as 



P* 



1 



l + 2 7 



67A* 1 - d 
1 + 2d 7 - q-^- 



a,* |2 (1 + 6 7 ) 2 



(3.21) 



l + 2 7 



67A* 
^| 2 ( 1 +67) S 



(3.22) 



P * xy (1 + 2 7 ) 2 



X* I-67 
<\2 (! + 6 7) 5 



(3.23) 
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FIG. 1: Shear-rate dependence of the elements of the reduced pressure tensor P*j = pj^/p for a = 1 in the three-dimensional 
case. The solid lines correspond to the results obtained for Model A (q = 0) while the dashed lines refer to the results derived 
for Model B for q — | by using the approximation (|3.16|l . 
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FIG. 2: Shear-rate dependence of the elements of the reduced pressure tensor P*j = Pff 1 lp for a = 0.5 in the three-dimensional 
case. The solid lines correspond to the results obtained for Model A (q = 0) while the dashed lines refer to the results derived 
for Model B for q = | by using the approximation (|3.16[1 . 

As expected, when the steady state condition (|2.4p applies locally, then d[ T = 0, A* = so that, according to Eqs. 
([3T7]) and ([3T8]l . the results for the pressure tensor are independent of the model interaction considered. Note that 
the same type of approximation (|3.16[) can be used to estimate the fourth-degree moments for q ^ 0, although this 
calculation will be omitted here for simplicity. 

Figures [T] and [2] show the dependence of the reduced pressure tensor PZ on the reduced shear rate a* for a 
three-dimensional system (d = 3) and two different values of the coefficient of restitution: a = 1 and a = 0.5. I have 
considered the exact results obtained for Model A and the results for Model B (when q — h)by using the approximation 
(I3.16p . It is quite apparent that, for a given interaction model, there is a weak influence of the inelasticity on the 
shear-rate dependence of the pressure tensor elements. Furthermore, for a given value of a, the differences between 
the results obtained for Model A and Model B are very small in the range of shear rates considered. This means 
that the rheological properties are quite insensitive to the interaction model considered, even for strong dissipation. 
The dependence of the pressure tensor on dissipation in the steady state [i.e., when the (reduced) shear rate and the 
coefficient of restitution a are not independent parameters but they are related through Eq. (|2.30|) ] is plotted in Fig. 
|3] along with the results obtained for IHS [23|, [27} ■ It can be observed that the IMM results reproduce very well the 
IHS predictions over the range of values of the coefficient of restitution analyzed. 
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a 

FIG. 3: Plot of the reduced elements of the pressure tensor P* as functions of the coefficient of restitution a in the steady USF 
state. The solid lines arc the results derived here for IMM while the dashed lines correspond to the results obtained for IHS in 
Ref. IH. 



B. First-order approximation 



The analysis to first order in the gradients is worked out in Appendix [Bl Only the final results are presented in 
this Section. The distribution function is of the form 



/ 



(i) 



X„ ■ Vn + X T • VT + X u : VSu, 



(3.24) 



where the vectors X„ and Xj- and the tensor X„ are functions of the true peculiar velocity c. They are the solutions 
of the following linear integral equations: 



X,, 



T 



2a 
dp 



(1 - nd n )P^ - C 



Xn 



(3.25) 



2 

dn 



) + (q+l)C + ac 



v dc a 



-C 



Xta = Yta, 



(3.26) 



dn 



aPf$+T()d T 



' dc x 



X u . 



a5kyX u ^xi — Y u . 



kt • 



(3.27) 



where Y„(c), Y-r(c), and Y„(c) are defined by Eqs. (|B10[) — ()B12|) . respectively. In addition, C is the linearized 
Boltzmann collision operator around the reference USF state: 



CX = - 



(j[/(°),X] + J[X, /(»)]). 



(3.28) 



It is worth noting that for q = |, Eqs. (|3.25p (|3.27|) have the same structure as that of the Boltzmann equation for 
IHS 12]. The only difference between both models lies in the explicit form of the linearized operator C. 
With the distribution function /W determined by Eq. (|3.24| . the first-order corrections to the fluxes are 



(!) _ _ dSu k 



dT dn 
lh, ''''Ik, 



(3.29) 
(3.30) 
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where 

Vijki = -I dcmciCjX Ui kt>(c), (3.31) 



f Tfi 

I dc —c 2 c i X T . j (c), 



(3.32) 



/*«=-/ dc—c 2 CiX, hj (c). (3.33) 

Upon writing Eqs. (|3.29p - (|3.33p use has been made of the symmetry properties of X n ^, Xt,i and X u jj. In general, 
the set of generalized transport coefficients rjijkt, K-iji and are nonlinear functions of the coefficient of restitution a 
and the reduced shear rate a*. The anisotropy induced in the system by the presence of shear flow gives rise to new 
transport coefficients, reflecting broken symmetry. The momentum flux is expressed in terms of a viscosity tensor 
"Hijkiip*!®) of rank 4 which is symmetric and traceless in ij due to the properties of the pressure tensor fy . The 
heat flux is expressed in terms of a thermal conductivity tensor (a* , a) and a new tensor /lij (a* , a) . Of course, for 
a* = and a = 1, the usual Navier-Stokes constitutive equations for ordinary gases are recovered and the transport 
coefficients become 

VijM — > m (fiikSje + 5 jk 5 ie - -^>ij$kt \ , — > KoSij, fiij — > 0, (3.34) 

where 770 — p/^o^and «o = d{d + 2)r)o/2(d — l)m are the shear viscosity and thermal conductivity coefficients for 
elastic collisions 



IV. SHEAR-RATE DEPENDENT TRANSPORT COEFFICIENTS 



This Section is devoted to the determination of the generalized transport coefficients rjijkt, Kij, and fiij associated 
with the momentum and heat fluxes. Let us consider each flux separately. 



A. Momentum flux 



To first order in the hydrodynamic gradients, the momentum flux is given by Eq. (|3.29|) . To get the coefficient 
rjijke, I multiply both sides of Eq. (|3.27|) by mciCj and integrate over c. The result is 

= -5 kl (1 - nd n ) Pjp - (s ik P$ + 5 jk P$) 

+ ll(P^-^ V M)d T P§\ (4.1) 
Upon writing Eq. (|4.1[) , use has been made of the results [HI, [I?} 

dc mciCjCXuM = -^o|2%w, (4.2) 



J dc mc lCj Y uM = -5 H (1 - nd n ) P$ ~ (s ik P^ + 5 jk P^ 

+ ^{P^-^yki)d T P^. (4.3) 

The generalized shear viscosity can be written as r\ij k i = VoVijkii -*) where r]ij k i( a *) is a dimensionless function of 
the reduced shear rate a* and the coefficient of restitution a. The dependence of T]*j k £ on temperature is through the 
reduced shear rate a* and so 

Td T r]ijke = Td T rior]* jke (a*) = (f - q)r)iju - qr]i]kea*da* h\if ijM (a*). (4.4) 
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Consequently, in dimensionless form, Eq. (|4.1|) yields 



\ a *Ky + C^j [1 - 9(1 + a* da-)) VijM - a * (StxVjykt + S JxV*yke ~ S kyV*jxi) ~ »0\2Vt 



-6 ke a*d a *P* - (S lk P* e + 8 jk P* t ) 
2 

"d 



+ {Pu ~ a'nlyu) (1 - W*d a ,)P*, (4.5) 



where use has been made of the identity 



d O (o) d , iS A „ <9 



In the absence of shear field (a* — 0), P*j — Sij, and so Eq. f|4.5|) has the solution 

^o|2 - C*(l - 9) f^fc^ + SjkSu - ^|<%£jw^ • (4-7) 

This expression coincides with the one derived by Santos [17J for IMM for vanishing shear rates. Beyond this limit 
case, in general Eq. (|4.5|) is a nonlinear differential equation that must be solved with the appropriate boundary 
conditions to get the hydrodynamic solution. The simplest model is that of Model A, in which case the set (|4.5j) 
becomes a set of coupled algebraic equations. 

1. Model A 

When q = 0, then the shear viscosity obeys the equation 



(l a * P xy ~ Vtjkt - a * ( S ixVjyki + S JxV*yke ~ SkyVijxi) 

= -8ua*d a *P* j - (S ik P* e + SjkPu) 

+~ d {Pu-<fr,t vU )P* r (4.8) 

The explicit dependence of r]*- M on a* and a can be obtained by solving the set of algebraic equations (|4.8[) . As an 
example, the coefficients of the form f]*j xy are obtained in detail in Appendix [C] The nonzero elements of Tj*j Xy are 
given by 

„• = — IZ*L (4 9) 

' xyxv w* |2 (l + 2 7 ) 2 (l + 6 7 )' [ 1 

„ 4 1 - d a 

Vxxxy ~^~ 2 d (l + 2 7 ) 3 (l + 6 7 )' ( ' j 



4 

cLJ* ]2 (1 + 27)8(1 + 67) 
When a = 1, Eq. (|4.9|) reduces to the one previously derived for elastic gases (29j . 



^yyxy Vzzxy j. .* , i " 

(4.11) 



2. Mode/ £ 

When g 7^ 0, the coefficients r/*j ke must be obtained numerically. However, assuming that the dependence of "q*- M 
on a* and a for q =/= does not differ significantly from the one obtained for Model A (q = 0), one can consider again 
the approximation 

1#( = V*jM,A + 1 A Vijke + ■■■ , (4-12) 
to get analytical results for r\*^ kl . In Eq. (|4.12[) . rj*- M A refers to the known result for Model A. The unknown A?/*^ 
can be determined by inserting (|4.12p into Eq. (|4. 5|) and retaining only linear terms in q. As an illustration, the 
quantities Ar]*j Xy are also evaluated in Appendix [C] Their explicit expressions are given by Eqs. (|C6[) - (|C10[) . 
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FIG. 4: Plot of the reduced element rf xyxy as function of the coefficient of restitution a in the steady USF state for a three- 
dimensional system. The solid line is the result derived here for Model A (g = 0) , the dashed line refers to the result for Model 
B for q = | by using the approximation (|4.12|l . and the dotted line corresponds to the result obtained for IHS in Ref. [13 ] . 



B. Heat Flux 



The heat flux is denned by Eq. (|3.30p in terms of the coefficients Kij , Eq. (|3.32|) , and /i^ , Eq. (|3.33|) . To get these 
generalized Navier-Stokes coefficients, let us introduce the tensors 



I dc —CiCjCkX Tt i(c), 



(4.13) 



/7JI 
dc —CiCjCkX nt i(c). 

The generalized thermal conductivity Kij and the new coefficient faj are given by 



(4.14) 



(4.15) 



Let us consider first the tensor Kijki- To get it, 1 multiply both sides of Eq. (|3 . 26[) by ^CiCjCk and integrate over 



velocity. After some algebra, one gets 



where [If 



"211 



+ (^(T&rPW) + (g + 1)C - i/ | 8 ) ««M 
-a (S 

ix^jkyi ~t~ &jx^iky£ ~t~ ^kx^ijyi) 

v 2\l — v 0\3 f c , s , t \ ' ' ' ■ I '\ ' 1 
d~+2 jk °ik K 3t + OijKkg) = --^OT-Nijkt 



(l + a)[5d + 4- (d+8)a] 



Sd 



-vq, ^o|3 = 2^012, 



(4.16) 



(4.17) 



N 



(0) 



-ijU ~ I dcc tCj c k c e f { >{c). 
Upon writing Eq. (|4.16p . use has been made of the results 



a m rv "2\l- ^o|3 

dc —aCjCkLXr/ = ~vo\z K ijki d + 2 — jk °ik K ji + OijKkt) , 



(4.18) 



(4.19) 
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FIG. 5: Plot of the reduced element r\y yxy as function of the coefficient of restitution a in the steady USF state for a three- 
dimensional system. The solid line is the result derived here for Model A (g = 0) , the dashed line refers to the result for Model 
B for q — | by using the approximation (|4.12|l . and the dotted line corresponds to the results obtained for IHS in Ref. 



m 

dc —c 2 CjC k Y T ,, 
2 



rn 
~2 



2mn 



(o) 



P^d T P, 



P (0)q p(0) 



(4.20) 



As happens in the case of the shear viscosity, the thermal conductivity tensor K.y can also be written as Ky = 
KoK*Ja,a*), where Kt^j is Si dimensionless function of a* and a. Thus, 

Td T Kij k £ = Td T K K* jke (a*) = (1 - q)n ijM - qK ljM a*d a * ln/t?- M (a*). (4-21) 
In addition, the derivative 9tN^^ of the fourth-degree moments of the zeroth-order distribution can be evaluated as 



Td T Nj 0) 



ijkt ~ m 2 

where I have introduced the reduced fourth-degree moments 

N *jke( a *) = 1 



nT 2 

4^(2-qa*d a *)N* m (a*), 



N 



(0) 



Finally, in dimensionless form, Eq. (|4.16[) becomes 



-a*P* + C [1 - 9(1 + a*d a ,)} k* + 



4 nT 2 ij W 



2a 



— (1 - ga*0 o .)i* + (9 + 1)C* - ^o|3 



2|1 ^0|3 



d + 2 
4(d-l) 



{2-qa*d a *)N* jH 
[P^(l-qa*d a ,)P* ( 



d(d + 2) 



(4.22) 



(4.23) 



(4.24) 
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FIG. 6: Shear-rate dependence of the ratio rj X y Xy {a*) / r)% yxy {p) for two values of the coefficient of restitution a in the three- 
dimensional case. The solid lines correspond to the results obtained for Model A (q — 0), while the dashed lines refer to the 
results derived for Model B for q — i by using the approximation (I4.12[l . 
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FIG. 7: Shear-rate dependence of the ratio Vyyxy( a *) /VyyxyiO) for two values of the coefficient of restitution a in the three- 
dimensional case. The solid lines correspond to the results obtained for Model A (q = 0), while the dashed lines refer to the 
results derived for Model B for q = i by using the approximation (|4.12[l . 



where u*^ 3 = v \ 3 /v and v*^ = u 2 \i/v . 

The set of equations defining the elements of the tensor /iijke can be easily obtained by following similar mathe- 
matical steps as those just made for the thermal conductivity tensor. Thus, the reduced tensor fJ-*j kf = (n/TKo)fJ-ijki 



is determined from 
'2 



) [2 - q(l + a*d a .)] H* M - {^{d a 'P*y) - C*) < jU 



2|1 0|3 / j- *ijr * i x * ^ 



d + 2 

4(d-l) 
d(d + 2) 

d-l 
d{d + 2) 

+P* k (l - a*d a *)P*t + f£(l - o*5..)Pjm] 



(l-a*d a *)N* jM 



(4.25) 

In the absence of shear rate (a* = 0), the solutions to Eqs. (|4.24|) and (|4.25|) give the explicit forms for the thermal 
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FIG. 8: Plot of the reduced elements n yy and K* y of the thermal conductivity tensor as functions of the coefficient of restitution 
a in the steady USF state for a three-dimensional system. The solid lines correspond to IMM for Model A (q = 0), while the 
dashed lines refer to the results for IHS derived in Ref. 



conductivity tensor K*j and the tensor ju?.. They can be written as 

4 = ^^jrC 1 + c ) (^11 - 2C*) _1 , (4.26) 



where the fourth- cumulant c is 111 



1 + cv, 



2|1 



(2 - q)C 



(4.27) 



c(a) 



12(1 - a) 2 



Ad -7 + 3a(2 - a) ' 



(4.28) 



Equations (|4.26[) and (|4.2T|) coincide with the previous expressions derived by Santos [l7j for the Navier-Stokes 
transport coefficients of IMM associated with the heat flux. 

As in the case of the shear viscosity, Eqs. (|4.24|) and (|4.25|) become algebraic for Model A (q = 0). Even for this 
model, although the solution to Eqs. (|4.24p and (|4.25|) is simple, it involves a quite tedious algebra due to the presence 
of the fourth-degree moments N*^ k[ whose expressions for hard spheres (d — 3) are displayed in Appendix [X] As an 
illustration, the detailed form of the elements n* y and n* y of the thermal conductivity tensor is obtained in Appendix 
[Cl A more complete list of the coefficients k*. and /x*- can be obtained from the author upon request. 



V. COMPARISON WITH THE TRANSPORT COEFFICIENTS FOR IHS 



As said in the Introduction, the expressions for the generalized transport coefficients of IHS described by the 
Bolzmann equation have been recently derived [l2|,[2(| when the steady state conditions (|2.30[) apply. These expressions 
have been obtained by using a BGK-like kinetic model of the Boltzmann equation for a three-dimensional system [l2l ]. 
In this Section, some of the coefficients obtained here for IMM will be compared with those presented for IHS in the 
steady state. Beyond the steady state conditions, the dependence of the generalized transport coefficients r]*j ke and 
K*j on both the reduced shear rate a* and the coefficient of restitution a will be also analyzed for several situations. 
In this case and to the best of my knowledge, there are no available IHS results and so a comparison between IMM 
and IHS cannot be carried out. 

Let us study first the coefficients Vxyyy an d Vyyyy m the steady USF state. As said before, in this situation a* 
is related to a through Eq. (|2.30[) and so, the coefficient of restitution is the only control parameter of the system. 
Figures S] and [5] show the a-dependence of vt vvv and Vyyyyi respectively, as given by Model A, Model B with q = | by 
using the approximation (|4.12[) , and IHS L 12j in the steady USF state. It is apparent that the agreement between the 
predictions of Model B for IMM (which tries to mimic the true IHS model) and IHS is excellent in the whole range of 
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FIG. 9: Shear-rate dependence of the ratios K, yy (a*) / n yy {Q) (solid lines) and K* y (a*)/n^. y (0) (dashed lines) for two values of 
the coefficient of restitution a in the three-dimensional case. All the results have been derived from Model A (q = 0). 

values of a analyzed. This confirms the reliability of IMM to reproduce the main trends observed for the coefficients 
Vijkf Regarding the results for Models A and B, we observe quantitative differences between both models for the 
coefficient Vyyyy f° r strong dissipation. In the case of general perturbations, the shear-rate dependence of the ratios 
Vy V xy( a *)/ r lyyxy(Q) an d r lxyxy{ a *) /VxyxyW is plotted in Figs. and [3 respectively, for two values of the coefficient 
of restitution: a = 1 and a = 0.5. Here, f]*j yy {0) is the value of T]*j yy for vanishing shear rates given by Eqs. (|C3|) 
for Model A and Eqs. (|C11[) and (IC12[) for Model B. We observe that, in general, the ratios r]i 3 xy( a *)/ r 1ijxy(Q) are 
monotonically decreasing functions of the shear rate (shear-thinning effect), except in the region of small shear rates 
in the case of Model B for Vx V xy( a *) / VxyxyW when a = 0.5. Figures [5] and [7] also show that, at a given value of a*, 
the value of T)*j xy (a*) /rj*j (0) decrease with dissipation. This means that the inelasticity produces an inhibition of 

the momentum transport since the value of for inelastic collisions is smaller than the one obtained in the elastic 

case. 

Next, the thermal conductivity tensor k*^ is considered when the steady conditions (|2.30p hold. The dependence 
of the elements k* y and k* on a is plotted in Fig [5] for IMM (in the case of Model A) and for IHS [e| • As happens 
for the usual thermal conductivity coefficient the trends observed for IHS are strongly exaggerated by IMM, 
especially at high inelasticity. It is possible that the disagreement between both interaction potentials at the level 
of the heat flux might be mitigated in part if one considered Model B instead of Model A for IMM. However, this 
would require to estimate the fourth-degree moments of the reference state by using the approximation (|3.16[) for 
5 = 5, which is quite an intricate problem. The shear-rate dependence of K*j is illustrated in Fig. [5] for the ratios 
K xy( a *) / K xy(®) an d K l v ( a *) I K yy(fy f° r two values of the coefficient of restitution: a = 0.5 and a = 1. Here, k*(0) 
is a Burnett coefficient given by Eq. (|C28|) and K* y (0) is the Navier-Stokes thermal conductivity coefficient given by 
Eq. (|4.26p [or, equivalently, Eq. (|C29|) ]. We see that n* y increases with a* in the region of shear rates considered 
for elastic collisions, while it does not present a monotonic dependence on a* in the inelastic case since it reaches 
a maximum and then decreases with respect to its Navier-Stokes value. Consequently, while for elastic gases the 
shear flow enhances the transport of energy along the direction of the gradient of the flow velocity (y axis), this is 
not the case for inelastic collisions for large shear rates, where the heat flux is inhibited by the shearing motion. The 
off-diagonal element k* measures cross effects in the thermal conduction since it gives the transport of energy along 
the x axis due to a thermal gradient parallel to the y axis. This cross coupling does not appear in the linear regime 
as the heat flux qjp is at least of Burnett order (proportional to a*dT/dy). The element n* xy is negative and its 
absolute value decreases with respect to its Burnett value K xy (0) as the shear rate increases, regardless of the value 
of dissipation. 

VI. DISCUSSION 

In this paper, the transport properties of <i-dimensional IMM in inhomogeneous states close to the USF state have 
been analyzed. The physical situation is such that the granular gas is in a state that deviates from the USF by small 
spatial gradients. Given that the system is subjected to a strong shear flow, the corresponding transport coefficients 
associated with the irreversible heat and momentum fluxes are nonlinear functions of the shear rate. The explicit 
evaluation of these coefficients has been the primary objective of this work. The search for such expressions has been 




18 



prompted by recent results [12j | derived from a simple kinetic model of the inelastic Boltzmann equation, which provides 
explicit expressions for the above coefficients in the case of IHS. Since this kinetic model is a simplified version of the 
nonlinear (inelastic) Boltzmann equation, a still open problem is to get the explicit shear rate dependence of these 
transport coefficients when the true Boltzmann collision operator is considered. For this reason, due to the technical 
difficulties associated with the mathematical structure of the Boltzmann equation for IHS, the problem studied in Ref. 
[l2l ] has been revisited here for IMM. For this interaction model, the collision rate becomes independent of the relative 
velocity of the two colliding particles so that the velocity moments of order k of the Boltzmann collision operator can 
be exactly written in terms of moments of order k' < k [l5| . Thanks to this feature, the second- and fourth-degree 
velocity moments corresponding to the (pure) USF state has been recently obtained . The knowledge of the latter 
moments provides the proper basis to explicitly evaluate transport properties around USF. 

As said before, I have been interested in a situation where weak spatial gradients of density, velocity, and temperature 
coexist with a strong shear rate. Under these conditions, the Boltzmann equation is solved by means of an extension 
of the Chapman-Enskog method to arbitrary reference states [2(J. In the case of the USF state, due to the anisotropy 
induced by the shear field, tensorial quantities are required to describe the momentum and heat fluxes instead of the 
usual Navier-Stokes transport coefficients [3,[T3]- In the first order of the expansion the momentum and heat fluxes are 
given by Eqs. (|3.29p and (|3.30p . respectively, where the components of the set of generalized transport coefficient i)ijhl, 
Kij, and are the solutions of Eqs. (|4.5p . (|4.24[) . and (|4.25p , respectively. Such coefficients are nonlinear functions 
of both the shear rate and the coefficient of restitution. As expected, there are many new transport coefficients 
in comparison to the case of states near equilibrium or near the cooling state. These coefficients provide all the 
information on the physical mechanisms involved in the transport of momentum and energy under shear flow. Of 
course, the usual form of the Navier-Stokes coefficients for IMM is recovered in the absence of shear flow [l7| • 

The purpose of this work has been two-fold. First, the evaluation of the shear-rate dependent transport coefficients 
of IMM is worthwhile studying by itself as a tractable model to gain some insight into the combined effect of shear 
flow and dissipation on the transport properties of the system. Second, the comparison between the exact results 
derived here for IMM with those previously obtained for IHS by using a simple kinetic model allows us to assess the 
relevance of IMM to reproduce the trends observed in the case of IHS. Recent comparisons [13, EH [HI carried out 
between both interaction models show good agreement, especially in the USF problem [j^ . 

The results derived here are given in terms of an effective collision frequency u>, which depends on space and time 
through its dependence on the density and temperature. Here, I have considered the general dependence u> oc nT q , 
where q = is closer to the original Maxwell molecules for elastic collisions and q = \ is closer to hard spheres. To 
make some contact with the results obtained for IHS [l2[ , w is chosen to reproduce the cooling rate ( of IHS evaluated 
at the local equilibrium approximation. This yields the relation (|2.4[) . While Model A (q — 0) lends itself to exact 
analytical results for rjijke, K ij, an d /t*y, that is not the case for Model B (q ^ 0) since the latter requires to be solved 
numerically to determine the dependence of the above generalized coefficients on both the shear rate and dissipation. 
However, some analytical approximate expressions can be obtained by considering an expansion of the transport 
properties in powers of the interaction parameter q. In this case, the results obtained for q = \ have been compared 
with those obtained for IHS. Figures [4) El [3 and [9] illustrate this comparison when the steady state conditions (|2.17|) 
apply locally [and so, a* and a are related through Eq. (|2.30|) ]. It is apparent that the IMM predictions compare quite 
well with IHS in the case of the viscosity tensor T)ijkt, although the discrepancies between both interaction models 
increase in the case of the thermal conductivity tensor Kij . In this latter case, the IMM results tend to exaggerate 
the trends observed for IHS. 

An application of the results obtained in this paper would be to perform a linear stability analysis of the hydro- 
dynamic equations with respect to the steady simple shear flow state. This analysis allows one to determine the 
hydrodynamic modes for states near USF as well as the conditions for instabilities at long wavelengths. Previous 
results for IHS [12j in the three-dimensional case indicate that the USF is linearly stable when the perturbations 
are along the velocity gradient (y direction) only, while it becomes unstable when the perturbations are along the 
vorticity direction (z direction) only. The question now is whether these conclusions are similar in the case of IMM. 
Another possible direction of study is to extend the analysis made here for a single gas to the important subject of 
granular mixtures. Previous works carried out by the author and co-workers (l8l.[l3| for mixtures of IMM have shown 
the tractability of the Maxwell kinetic theory for these complex systems and stimulate the performance of this study 
in the near future. 
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APPENDIX A: FOURTH-DEGREE MOMENTS IN THE USF STATE FOR MODEL A 



In this Appendix, the explicit expressions of the set of the fourth-degree moments of the zeroth-order distribution 
f(°' are given in the case of Model A (q = 0) for a three-dimensional system (d — 3). As the set of independent 
moments, it is convenient to take the (dimensionless) Ikenberry moments [3lj 



where 



{M^ Q , M^^, M^^, MQ^yy, ^^Q^ZZZ^ ^^^XfJ^ ^Q^^y, ^0 \ X J/ J/ J/ } ) 



ic 2 J ) /(»I(c), 



dc c ( CjCj 



(Al) 



(A2) 



(A3) 



M, 



o\ijke 



dc 



1 / my 
n V2T/ 



CiCjCkCi - -c 2 (ciCjSki + CiCkSjt + CiCeSjk + CjC^Su + c ci8 ik + CkCe,8 i: j) 



/ (0) (c). 



By using matrix form, the moments (|Alj) are given by [161 ] 

m = c 1 e, 

where Al is the column matrix defined by the set (|A1|) 



M 



M. 



2\xx 



M 2\yy 

M* 
0\yyyy 

M*. 

lvl 0\zzzz 

M* 



2\xy 
)x 

\ 1V1 0\xyyy 



M* 

0\xxxy 



and C is the square matrix 



where I is the 8x8 identity matrix and 



C = 4w*| 27 X + £', 



£' 



( 


W 4|0 














4a* 








\ 







W 2|2 











32 * 
21 " 


2a* 















W 2|2 








21 " 





2a* 















W 0|4 





150 * 
245" 





-fa* 


















W 0|4 


^-a* 

245" 


fa* 


fa* 






15" 


fa* 


fa* 


-2 a * 


-±a* 


W 2|2 















±^a* 

49 " 


49 " 


-§<z* 


"Tl"* 





W 0|4 







V 





49 " 


±^a* 

49 " 


2a* 


fa* 








<\4 


J 



with [If 



(l + a) 2 (5 + 6a -3a 2 ) 



J 4\0 



48 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 
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The column matrix C is 



where 



W 2|2 - 



(l + a) 2 (34 + 21a-6a 2 ) 
168 : 



_ (1 + a) 2 (150 + 21a-3a 2 
o|4- 



/ Ci \ 

c 2 
c 4 

Co 



d = 1\\ - \; (6M* [yy 2 + 2M * |o 



2 A 3 M 0| TO - 3 A 4 ( M 0\xy 2 ~ 3M 0|yy 2 ) . 



C4 35 A * (' 



81M *| W 2 - 48M %„ 2 



Cs - 35^ ( 



I S1M* |W 2 + 12M % tf a 



^6 — 2 A 3"^oW + ^i^0\xy^0\w' 



36 

d = —^-\M Q \ xy M ^ yy , 



7 /v 5 J "0|a;y- IM 0|yjf 



C 8 = — AjM n *u..Af, 



In Eqs. (|A13|) — (|A20[) . I have introduced the second-degree moments 



where P* is given by Eqs. (|2.22j) - (|2.23p and the quantities 



A * = 144 (1 + a)2(ll_6a + 3a2) 



\* — 

A 2 — 



(l + a) 2 (l + 6a-3a 2 ) 



(A10) 
(All) 



(A12) 
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(A13) 
(A14) 
(A15) 
(A16) 
(A17) 
(A18) 
(A19) 
(A20) 

(A21) 

(A22) 
(A23) 
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A;= (1 + a)2(22 - 21 " + 6 " 2) . (A24, 



„ (l + a) 2 (21a-l-6a 2 ) fl 
K = LA ^ l L , (A25) 



(1 + a) 2 (39 - 21a + 3a 2 ^ 
378 



A 5 = " -■ (A26) 



The explicit dependence of the fourth-degree moments (|A1|) on a* and a can be determined from Eq. (|A5[) . These 
expressions are very long and not very illuminating so that they will be omitted here for the sake of clarity. 

The knowledge of the fourth-degree moments of the reference USF state is needed to compute the thermal conduc- 
tivity tensor and the tensor ^ [see for instance, Eqs. (|4.24p and (14.25[) ]. The relationship between the (canonical) 
moments defined in (|A1[) and the moments N*- M defined by Eqs. (|4.18|) and (|4.23p is given by 

Kxxx = \m* Mq + ®M* [xx + -{M^ yym + M* ]zzzz ), (A27) 

Kvvv = \ M h + 6 j M 2\vv + M oW (A28) 
NLyy = ^M: l0 + ~ (M* ]xx + M; lyy ) - 7 -M* [yyyy - \M^ ZZZZ) (A29) 
1 1 1 * 7 * 

N XX zz = 15^4|0 _ j^2\yy ~ Q^Q\yyyy ~ qMq\zzzz> (A30) 



Nyyzz — — ^4|0 _ j^2\xx + Q^0\yyyy + qMq\ zzzz , (A31) 



3 

^xxxy = j^2\xy + ^0\xxxy^ (A32) 

3 

Nxyyy = y A^2|a;y + Af |a:yj,y ; (A33) 

A£y„ = jM* {xy - (M* {xxxy + M* ]xyyy ). (A34) 

APPENDIX B: CHAPMAN-ENSKOG-LIKE EXPANSION 

In this Appendix, some technical details of the Chapman-Enskog-like type of expansion used in the main text are 
described. First, the action of the operators <9 t on the hydrodynamic fields can be obtained from the corresponding 
balance equations associated with the Boltzmann equation (|3.2p . They are given by 

d t n + u Q ■ Vn = -V • (n5u), (Bl) 
d t Su + a • Su + (u + Su) • V<5u = -(mnj^V • P, (B2) 



i ndtT + ^ n (u + Su) • VT + aP xy + V • q + P : VSu = ~pC, ( Li:', ) 
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where the pressure tensor P, the heat flux q and the cooling rate £ are defined by Eqs. (]2.11[) — (|2 . 13[) . respectively, 
with the replacement V — > c. From Eqs. (|Bl[) - (|B3p . it is easy to see that the action of the operator is given by 
Eqs. (|3.9|) and (|3.10p . while in the case of the operator <9 t the result is 

d t (1 W u • Vn = -V • (n5u), (B4) 

a (1) (5u+ (u + Su) ■ V5u = — — V • P (0) , (B5) 

mn 

ind^T + ^n{u + Su) ■ VT = -aP£) - P (0) : VSu, (B6) 
2 2 " 

where 

Pjp = f dcmciCjf^ic). (B7) 
Inserting the expansions (|3. 5j) and (|3.6|) into Eq. (|3.2p . one gets the kinetic equation for the velocity distribution 

d 



d * ~ aVy w x +c ' = " dt + (v + Uo) ■ v /(0) ' (B8) 



where C is the linearized Boltzmann collision operator around the USF defined by Eq. (|3.28[) . Use of Eqs. (|B4p ~(|B6 
in Eq. ||B8j) yields 



( 5 * 0) ~ aVy W + C ) = Y » ' Vn + Yt • VT + Y„ : V<Su, 



(B9) 



where 



5/(0) IdfMdpW 

'~a — Ci + ~iTx « > ( B1 °) 

On p oouj on 



= n ^T^ ~ ^ + ( p « " a M • (B12) 

Equations (|B9)) - (|BT2)) are similar to those obtained for IHS [II], except that C (1) = for IMM. The solution to Eq. 
(IB9I) is of the form given by Eq. (|3.24[) . To get the corresponding integral equations verifying the unknowns X n ,i, 



Xt,i, and X u ^j, one has to take into account the action of the time derivative S t on the temperature and density 
gradients, 

df ] V t T = V;«9 t (0) T 

= f^(l " "Wi? —) V ^ ~ ( + (« + 1)0 V»T, (B13) 
* dn z * n J \ an 1 



d[ 0) W t 5u = Vid^Suj = -a jk Vi5u k . (B14) 

Substituting the expression of given by (|3.24| into ()B8() . the integral equations for X n) j, -X^, and X U) ij are 
identified as the coefficients of the independent gradients. This leads to Eqs. (|3.25p - (|3.27p . 
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APPENDIX C: TRANSPORT COEFFICIENTS 



In this Appendix, some of the reduced transport coefficients rj*j ki , k*j, and \i*j are explicitly computed. Let us 
consider first, the coefficients of the form f]\^ y for Model A. According to Eqs. (|4.8| . these coefficients verify the 
equations 



^ P xy H)|2 ) Vxy 



a ^yyxy 



o 

_ p*2 p* 

^ xy yy> 



(CI) 



f— P* -LO*,\n* +— P* 17* = -P* P* 

\ d y ' 2 / yyxy d w lx v x v j w 



(C2) 



where the expressions of the reduced elements P?- are given by Eqs. (|2.22p and (|2.23p . Since JlL^ = Vzzxy = ■■■ = 
Vddxy an<1 ^ijW * s a t race l ess tensor, then rf xxxy = — (d — l)r]* yxy . The solution to the set of equations (|C1[) and 
(|C2j) leads to Eqs. (|4.9p and (|4.1ip when one takes into account the relations (|2.25p and (|2.26|) . For small shear rates 
(a* -c 1), the coefficients Vijxy behave as 



Vxyxy 



12, 



4 1 - 5j x d^ 



nxy .,* 



(C3) 



For Model B, one can estimate T]*j by using the approximation (|4.12|) . Neglecting nonlinear terms in q, it is easy 
to see that the quantities Ar]*j xy are determined from the set of coupled equations 



4a* 



—f P xy,A ~ w 0|2 ) A V*xyxy ~ a * A V*y 



j \ AP xy P xy,A a Vxyxy ,A AP xy) 



AP* 



+ ( a * P xy,A + C) (1 + a*9a*)Vxyxy,A 



J^xy,A a Vxyxy.AJ a ^a*P X y,A^ 



(C4) 



2n* 

—j-^xy,A - ^ 



012 



A n 



y <)■•'!> 



In* 2 

p* A„* - —AP* ( P* — n*n* \ 

^ W,A Ixyxy W w ^ 'lxyxy,AJ 

+ (a*Pxy,A + C) (1 + O^.O^.A 



j { P xy,A a Vxyxy,Aj a ^ a * P yy, 



A- 



(C5) 



where the elements APy are given by Eqs. (|3.17|) and (|3.18|) and the viscosities r)* y A and T]* A f° r Model A 
are given by Eqs. (|4.9[) and (|4.1ip . respectively. The solution to Eqs. (|C4[) and (|C5p . along with the constraint 
Arj* xxy = — (d— l)Ar]* yxy , provides the explicit form for the terms Arf-j xy . Their nonzero elements can be written as 



A V[ 



2K 



""" da;*2 2 (l + 27) 2 (l + 67) 4 ' 
22(1 - d)A yy 



d 2 w*f 2 (l + 2 7 )3(l + 67) 4 ' 



(C6) 
(C7) 



A Vyyxy A Vzzxy 



2aA, 



d 2 w*f 2 (l + 2 7 )3(l + 67) 4 



where 



A xy = 72d^, 2 7 4 - 12^ 9dcJo| 2 + 2(5 + 2d)( 



■2^ 



2(13 - 2d)C - 17dw*, 2 



(C8) 



3dwo| 2 + 2(ll + 4(i)C* 



(C9) 
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A yy = 144dw* |2 7 3 - 12 7 2 2dwo| 2 + (14 + 5d)C 
+ (d+10)C*. 

For small shear rates, the coefficients Vijxy have the expansions 



47 



6dw5 ]a + (4 + 7d)C 



txyxy 



'02 



1 - 2q 



du> 



02 



02 



2q 



2(25 + 4d)C* + 3dw5 
dc^ 2 



12 



-(1-S jx d) [1-qC 



10 



2dcJn 



(CIO) 

(Cll) 
(C12) 



The evaluation of the coefficients k\- and /x*^ is more intricate than that of 7]*j ki since it involves the fourth-degree 
moments of the zeroth-order distribution. Here, only the coefficients k* and k*, for Model A in the case of a 

j yy xy 

three-dimensional system (d — 3) will be considered. These coefficients are given by 



t-* _ ^.* _[_ If* I ty* 

xy xxxy ' yyxy ' zzxyi 



(C13) 
(C14) 



In the case of Model A (q = 0), the coefficients K* jke obey the set of coupled equations 



^*Ky+C 



y a\3 



U 2\l ^013 



a * {pixK*jkyl + 3jx K ikyl + ^kx K ijyl) 

2 / 

/ Q ?\r* P* P* P* P* P* P* \ 

"Yg \° Iy ijke ~ r kj r il ~ r ik r ji ~ r ij r kl) ■ 



(C15) 



Equation (IC15[) shows that K* y and K* y are coupled, so that their calculation involves the set of coefficients 

K^xxyyi ^yyyy ^zzyyi ^xyyy ^xxxyi ^yyxyi ^zzxy\' (^1^) 

By using matrix notation, the set of seven coupled algebraic equations for the unknowns (|C16j) can be written as 



Vaa'Qa'=Tla, <T = !,■■■, 7. 

Here, Q is the column matrix defined by the set (|C16[) . "P is the square matrix given by 



(C17) 
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1 + a 


[127(1 + a) 






\ a * P xy + 
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" ^0 3 


12 


+ 7a 


f 



(C18) 



The column matrix is 



x 0|3 v 1\\ 

5 

/ 87V. 



(C19) 
(C20) 
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15 



o p*2 

^ r xy 



P* P* \ 

xxyy xy xx yy 

Q 7V"* o p*2 

yyyy yy 

Q AT* p*2 

° yyzz zz 
q TV"* Q p* p* 

xyyy yy xy 

Q TV* Q p* p* 

c xxxy ° xx xy 
Q TV* Q p* p* 

xyyy ° £ yy 1 xy 

C/V"* p* p* 

xyzz zz xy 



(C21) 



25 



The solution to Eq. (|C17j) is 

Qa = {V- 1 ) aa .Tl al . (C22) 

This relation provides an explicit expression for the coefficients of the form K*j yy and rej,-™- From these expressions 
one can get the elements K* y and n* y . They can be written as 

K xy — ^ J K yy ~ ^ > (ozoj 

where 

= a*/3(3/3 + 80)^1 + ^-^2 [2a* 2 (3/3 + 20) + /3 2 (/3 - 20)] + ^ft 3 [/3 2 (/3 - 20) - 12a* 2 0] 

+2K 4 a* 2 (3/3 + 80) + (TZ 5 +K 6 + TZ 7 ) [/3 2 (/3 + 50) + 6a* 2 0] , (C24) 

# w = /3 2 (/3 + 50)^1 + ^ 2 [2a* 2 (/3 + 70) + /3 2 (/3 + 50)] + TZ 3 [(3 2 ((3 + 50) + 4a* 2 0] 

+27e 4 a*/3 (/3 + 50) - 2 (ft 5 +TZ 6 + TZ 7 ) a*/3<i>, (C25) 

A = 2a* 2 ( 6/3 + 230) + /3 2 (/3 2 + 10/30 + 2 5 2 ) . (C26) 
In the limit of small shear rates (a* — » 0), the coefficients K* y and n yy behave as 

< y « + <?(a* 3 ), « 4» («) + °(«* 2 )> ( C27 ) 

where 

(1) _ 112 162a 7 - 585a 6 + 981a 5 - 111136a 4 + 7272a 3 + 31267a 2 - 49855a - 9466 
Kxv {a > ~ 5 (1 - 9a) 2 (l + a) 3 [3a(a - 2) - 5] [3a(2a - 7) - 34] ' ( ' 

K (0) (a)= _J^ 17 + 9a(a-2) 

^ 1 > ~ 1 + a (9a - 1) [3a(a - 2) - 5] ' [ ' 

The coefficient K^ y (a) is a Burnett coefficient while K^/(a) gives the Navier-Stokes thermal conductivity coefficient 
for a three-dimensional inelastic gas, Eq. (|4.26[) . When a = 1, one recovers well-known results [4J for elastic collisions, 
namely, K y ° y = 1 and K^y = — \- 



[1] Goldshtein A and Shapiro M 1995 J. Fluid Mech. 282 75 

Brey J J, Dufty J W and Santos A 1997 J. Stat. Phys. 87 1051 

van Noije T P C, Ernst M H and Brito R 1998 Physica A 251 266 
[2] Brilliantov N and Poschel T 2004 Kinetic Theory of Granular Gases (Oxford: Oxford University Press) 
[3] Brey J J, Dufty J W, Kim C S and Santos A 1998 Phys. Rev. E 58 4638 

Brey J J and Cubero D 2001 Granular Gases Poschel T and Luding S (ed) (Berlin: Springer Verlag), Lectures Notes in 

Physics 564, pp 59-78 

Garzo V and Montanero J M 2002 Physica A 313 336 
Garzo V and Dufty J W 2002 Phys. Fluids 14 1476 
[4] Chapman S and Cowling T G 1970 The Mathematical Theory of Nonuniform Gases (Cambridge: Cambridge University 
Press) 

[5] Brey J J, Ruiz-Montero M J and Cubero D 1999 Europhys. Lett. 48 359 
Garzo V and Montanero J M 2002 Physica A 313 336 
Montanero J M and Garzo V 2003 Phys. Rev. E 67 021308 

Montanero J M, Santos A and Garzo V 2005 Rarefied Gas Dynamics 24 762 Capitelli M (ed) (American Institute of 

Physics) pp 797-802 and arXiv: |cond-mat/0411219| 

Brey J J and Ruiz-Montero M J 2004 Phys. Rev. E 70 051301 

Brey J J, Ruiz-Montero M J, Maynar P and Garcia de Soria M I 2005 J. Phys.: Condens. Matter 17 S2489 
[6] Cercignani C 1988 The Boltzmann equation and Rs Applications (New York: Springer Verlag) 



20 



[7] Garzo V and Santos A 2003 Kinetic Theory of Gases in Shear Flows. Nonlinear Transport (Dordrecht: Kluwer Academic) 
[8] Bird G A 1994 Molecular Gas Dynamics and the Direct Simulation Monte Carlo of Gas Flows (Oxford: Clarendon) 
[9] Brey J J, Dufty J W and Santos A 1999 J. Stat. Phys. 97 281 

[10] Tij M, Tahiri E E, Montanero J M, Garzo V, Santos A and Dufty J W 1999 J. Stat. Phys. 103 1035 

[11] Campbell C S 1990 Ann. Rev. Fluid Mech. 22 57 
Goldhirsch I 2003 Ann. Rev. Fluid Mech. 35 267 

[12] Garzo V 2006 Phys. Rev. E 73 021304 

[13] See for instance, Bobylev A V, Carrillo J A and Gamba I M 2000 J. Stat. Phys. 98 743 

Carrillo J A, Cercignani C and Gamba I M 2000 Phys. Rev. E 62 7700 

Ben-Naim E and Krapivsky P L 2000 Phys. Rev. E 61 R5 

Cercignani C 2001 J. Stat. Phys. 102 1407 

Ernst M H and Brito R 2002 Europhys. Lett. 58 182 

Ernst M H and Brito R 2002 J. Stat. Phys. 109 407 

Ernst M H and Brito R 2002 Phys. Rev. E 65 040301 

Baldasarri A, Marconi U M B and Puglisi A 2002 Europhys. Lett. 58 14 

Ben-Naim E and Krapivsky P L 2002 Phys. Rev. E 66 011309 

Krapivsky P L and Ben-Naim E 2002 J. Phys. A: Math. Gen. 35 L147 

Ben-Naim E and Krapivsky P L 2002 Eur. Phys. J. E 8 507 

Marconi U M B and Puglisi A 2002 Phys. Rev. E 65 051305; 66 011301 

Bobylev A V and Cercignani C 2003 J. Stat. Phys. 110: 333 

Bobylev A V, Cercignani C and Toscani G 2003 J. Stat. Phys. Ill 403 

Santos A and Ernst M H 2003 Phys. Rev. E 68 011305 

Bobylev A V and Gamba I M 2006 J. Stat. Phys. 124 497 

Ernst M H, Trizac E and Barrat A 2006 J. Stat. Phys. 124 549 

Ernst M H, Trizac E and Barrat A 2006 Europhys. Lett. 76 56 

Barrat A, Trizac E and Ernst M H 2007 J. Phys. A: Math. Theor. 40 4057 
[14] Kohlstedt K, Snezhko A, Sapozhni kov M V, Arn arson I S, Olafsen J S and Ben-Naim E 2005 Phys. Rev. Lett. 95 068001 
[15] Garzo V and Santos A 2007 arXiv: |cond-mat/0703704| 
[16] Santos A and Garzo V 2007 arXiv: cond-mat/0706.0475. 
[17] Santos A 2003 Physica A 321 442 
[18] Garzo V 2003 J. Stat. Phys. 112 657 
[19] Garzo V and Astillero A 2005 J. Stat. Phys. 118 935 
[20] Lutsko J F 2006 Phys. Rev. E 73 021302 

[21] Bobylev A V, Carrillo J A and Gamba I M 2000 J. Stat. Phys. 98 743 

Carrillo J A, Cercignani C and Gamba I M 2000 Phys. Rev. E 62 7700 

Bobylev A V and Cercignani C 2002 J. Stat. Phys. 106 547 
[22] Astillero A and Santos A 2007 Europhys. Lett. 78 24002 
[23] Santos A, Garzo V and Dufty J W 2004 Phys. Rev. E 69 061303 
[24] Lees A W and Edwards S F 1973 J. Phys. C 5 1921 

[25] Dufty J W, Santos A, Brey J J and Rodriguez R F 1986 Phys. Rev. A 33 459 
[26] Ikenberry E and Truesdell C 1956 J. Rat. Mech. Anal. 5 1 

Truesdell C 1956 J. Rat. Mech. Anal. 5 55 
[27] Brey J J, Ruiz-Montero M J and Moreno F 1997 Phys. Rev. E 55 2846 
[28] Santos A and Astillero A 2005 Phys. Rev. E 72 031308 

Astillero A and Santos A 2005 Phys. Rev. E 72 031309 
[29] Lee M and Dufty J W 1997 Phys. Rev. E 56 1733 
[30] Garzo V, 2007 J. Stat. Mech. P02012 

[31] Truesdell C and Muncaster R G 1980 Fundamentals of Maxwell's Kinetic Theory of a Simple Monatomic Gas (Academic 
Press: New York) 



